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The possibility of non-adiabatic electron pumping in the system of three coupled quantum dots 
attached to the leads is discussed. We have found out that periodical changing of energy level 
position in the middle quantum dot results in non zero mean tunneling current appeared due to 
non-adiabatic non-equilibrium processes. The same principle can be used for fabrication of a new 
class of semiconductor electronic devices based on non-stationary non-equilibrium currents. As an 
example we propose a nanometer quantum emitter with non-stationary inverse level occupation 
achieved by electron pumping. 
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I. INTRODUCTION 

Electron pumping in nanoscale structures attracts 
much attention nowadays [l],,2j],@,[3], [1], A great 
deal of the previous research works have been devoted to 
adiabatic electron pumping, the idea discussed by Thou- 
less [7] rather long time ago. The first (to our knowledge) 
experiment on electron pumping in single electron device 
was described in Then experiments in this direction 
were continued in a three-junction geometry by Pothier 
et al. [1] . Two phase shifted rf signals were used to real- 
ize a single electron pump: a device with current I — ef 
at zero bias voltage (/-frequency of the rf signal). Adi- 
abatic charge pumping based on periodical variation of 
the potential barriers forrned by the finger gates was also 
recently investigated in [l|. In these systems quantized 
current is connected with periodic adiabatic changing of 
the population of the quantum dot. 

Proposed in a number of papers photo-assisted tun- 
neling through coupled quantum dots @,[lo'|, [TT'l, fT2j 
is also an example of an electron pump. Pumping effect 
is achieved by applying an oscillating signal to the gate 
electrode or by irradiating the structure by monochro- 
matic [ll] and pulsed [11] microwaves. 

It was understood that for practical realization of 
quantized electron pump the phenomenon of Coulomb 
blockade is very important [l3| . General approach to the 
pumping through the interacting quantum dots in this 
regime is based on supposition that the charge relates 
to instantaneous chemical potential of a dot [15]. Using 
Coulomb blockade ideas another class of non-adiabatic 
quantized pumping of electrons in hybrid normal metal- 
superconductor structures was proposed in [l^, [l7| . 



These systems are more like a "turnstile" rather than a 
"pump" because quantized current directly proportional 
to gate frequency appears at finite (nonzero) value of ap- 
plied bias. But at low temperatures these systems are 
very promising as current standards [isj . 

Adiabatic charge pumping through three tunnel- 
coupled quantum dots attached to electron leads in the 
regime of strong Coulomb blockade was investigated in 
[19| . Slow variations of coupling strength between the 
dots lead to adiabatic changes of energy levels in the sys- 
tem. Time dependent charge redistribution caused by 
energy levels changes results in non-stationary adiabatic 
tunneling current. 

In the present paper we suggest a new type of electron 
pumping also in a system with three quantum dots, but 
based on non-equilibrium non-stationary tunneling cur- 
rents. The proposed device requires only a single ac gate 
signal contrary to other semiconductor devices which re- 
quire at least two rf signals with a definite phase shift. 
Mean current appears in our model due to non-adiabatic 
changing of electron level in a single quantum dot. 



II. THREE DOT MODEL OF ELECTRON PUMP 

We investigate non-stationary currents which flow in 
a three dot system shown in FigHJ The left and right 
dots have energy levels £2 and £3 constant in time. And 
the level position of the middle dot £1 is modulated by 
external gate voltage. 
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Quantum dots with energy levels £2 and £3 are also 
coupled to continuous spectrum states - massive leads. 
Hamiltonian of the system under investigation has the 
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FIG. 1: Schematic diagram of the three coupled quantum 
dots with energy level position in the middle quantum dot 
depending on the time. 



form: 



also presented in the equivalent differential form (except 
for the point t — t'): 
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Consequently, retarded Green's function which deter- 
mine spectrum re-normalization due to the tunneling be- 
tween the quantum dots can be written in the following 
form: 
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T12, T'13 are tunneling transfer amplitudes between 
the quantum dots and amplitudes Tk and Tp correspond 
to the tunneling processes between the quantum dots 
and continuous spectrum states, /a and c^(^p) / Ck{p)- 
electrons creation/annihilation operators in the quan- 
tum dots localized states and in the continuous spectrum 
states correspondingly. Energy values satisfy the follow- 
ing ratios: £2 > £_f and £3 < ep- Pumping effect appears 
if gate voltage switches the level £1 between energies £2 
and £3. 

In order to develop a theory for non-stationary cur- 
rent let us first describe in details charge relaxation pro- 
cesses in the system if we assume that at the initial mo- 
ment all charge density in the system is localized in the 
first (middle) quantum dot and has the value rii(O). At 
the first step we need to calculate exact retarded Green 
functions of the system. In the absence of tunneling be- 
tween the three quantum dots Green functions Gf^ {t~t ) , 
G^(t — t ) and G§^{t — t ) are equal to: 
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where 72 — t^v^T"^ and 73 — T^VpT^ are tunneling re- 
laxation rates from leftmost and rightmost dots respec- 
tively to the leads. Exact retarded electron Green's func- 
tion Gil i'^ the first quantum dot can be found from the 
integral equation: 



Ljii—(jii + i-rii Ji2'-^22'-^ll + "-^11 -'l3'-^33'-^ll \'^) 



Acting in 

"-^11 ''-^22 '"-^33 



turn by inverse operators 
this integral equation can be 



Gfi(i,i') ze(i -t')(Aie-*-^i +A2e-*^^(*^*') + 
+ A3e-*^^(*-*')) (5) 

Where eigenfrequencies £^1,2,3 can be found from equa- 
tion (see (g])): 

{E - £1) • (£; - £2 + 172) ■ {E- 83 + 173) - 
-Tf^ ■{E-e3 + 173) - Tfg • (£; - £2 + ^72) = 

(6) 

And coefficients Ai can be evaluated using integral 
equation for G{\: 
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where £» = £i - i-fi 

Further on we assume for simplicity that T12 = T13 = 
T. If £1 = £2 and £1 — £3 S> T, 7 coefficients Ai has the 
form: 
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El ~E2 ' 

If we disconnect the third quantum dot then coeffi- 
cients Ai and A2 give an exact solution for a system of 
two coupled quantum dots for all energy values £1 and 

£2- 
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Time evolution of electron density in the middle dot is 
determined by the Keldysh Green function [20j : 

G<,it,t) =ini{t) (9) 
Equations for the Green functions has the form: 

^^0-1 rp2 r<R 7^2 r<R \r<< — 



with coefficients F21, F31: 
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Now we can find all coefficients in : 
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(10) 

If G< (0,0) = mj^(e2) ^ 0, G<3(0,0) = mj^(£3) - 1 
and G^j(0,0) — ni(0) then Green function Gfi{t,t) is 
determined by the sum of homogeneous and inhomoge- 
neous solutions. Inhomogeneous solution of the equation 
can be written in the following way: 



G< (t,t') = iTfg [ dti [ dt2Gf^{t-ti) X 
Jo Jo 

xnf(£3)e^^^~^(*i-*^)Gi\(i2-t') 

(11) 

Homogeneous solution of the differential equation has 
the form: 



- iG<,it,t) = hit)e-''^^' + /2(t')e-'^^* + /3(i')e-^^^* 

(12) 

Function G^{t,t ) satisfies the symmetry relations: 
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Finally, time dependence of the filling number in the 
middle quantum dot ni(t) can be written as: 
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(13) 



[G<,{t,t)]* = -G<,{t,t) 
Then coefficients fi{t ) can be written as: 
flit) = Ae'^'^*' + Be'^'^'' +Xe*^3*' 
f2{t) = Ce'^'^*' + B*e'^'^*' + De'^'^' 
hit) = X*e*^i*' + D*e''^i^' + Ze*^3'*' 

(14) 

Since the solution has to satisfy homogeneous integro- 
differential equation, we are able to determine all coef- 
ficients. After some calculations we obtain that the fol- 
lowing proportionality takes place: 

f2it)^F2lhit) 

hit) = F^ihit) 



We see that there are six typical time scales in the 
considered system, which are described by the expression 
([T7| . Three of them we can identify as three relaxation 
modes with rates 2ImEi, 2ImE2 and 2ImE3 . Three other 
time scales are determined by the expressions Re(Ei — 
E2), Re(Ei - Eg) and Re(E2 - E3). These time scales are 
related with charge density oscillations between quantum 
dots, if the following ratio between Tij and 7 takes place: 
ry/7 > l/x/2. 

If we neglect for a moment the tunneling from the mid- 
dle dot to the right one, then for the system of two cou- 
pled quantum dots (ri3 = 0) only three time scales ap- 
pear and the equations (|15p . (IT6l) . pT)) are transformed in 
the following way: 
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Time dependence of the filling numbers in the first 
quantum dot ni(t) can be written as: 
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where coefficients A, B and G are determined as: 
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and coefHcients X, Z and D are equal to zero. 

In this situation we can distinguish two relaxation rates 
7res and ^nonres which characterises charge relaxation 
through an intermediate quantum dot in resonant and 
nonresonant cases: 



concrete system the accuracy of this approximation can 
be easily estimated from the exact equations. 

Pumping of electrons takes place if energy level ei{t) 
is a function of time and changes periodically (Fig. [IJ. 
We shall describe the most favorable case with T ^ 7. 
For current calculation we consider the situation of pe- 
riodically switching the position of level £1 by external 
(20) gate: 

ei(t) — £3 in the interval < i < to it means resonant 
tunneling between energy levels £1 and £3 

£i(i) — £2 in the interval < t < 2to - resonance 
between energy levels £1 and £2 

Time evolution of local electron density ni{t) in the 
central quantum dot can be determined from equation 

m (FigEJ. 
When 0<t<to 
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As we consider £1 — £3 > T, 7 then ^res 3> "inonres and 
small parameter ^nonres/lres exists in the theory. This 
allows us to write the following approximate relations for 
the system of three coupled quantum dots in the case 
£1 ~ £2 valid in the first order of the small parameter 
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and when Iq < t < 2to 
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Taking into account periodicity condition ni(2to) 
n^j one can find n\: 
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When £1 = £2 and £1 — £3 ^ T, 7 the exact equations 
T5|) can be transformed in the following way . 
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So, coefficients D, Z and X which are resposible for 
the "reverse" current to the right lead are much smaller 
than A, B and C, which correspond to "direct" current to 
the left, due to the appearance of the parameter -fj^zr^^- 
For simplicity we omit the terms with coefficients D, Z 
and X in equation (jl7p which determine time evolution 
of localized charge in the middle quantum dot. For any 



Results for ni{t) are shown in Figl2] Situation when 
frequency ft = l/2to of the level £i(t) switching is higher 
than tunneling rates 7res,7 is depicted by grey line. 
Black line corresponds to the case when frequency 51 is 
lower than 7resj7- It is clear that with the increasing 
of frequency the value of ni{t) always tends to the value 
1/2 and is almost independent on time at high gate fre- 
quencies. 

For low frequencies if presto ^ 1 the value is al- 
most equal to 1. The energy level £1 is filled up almost 
completely during the pumping cycle (for considered sit- 
uation when energy level £2 is well above and energy 
level £3 is well below the Fermi energy). Non-stationary 
tunneling current through the system appears for zero 
applied bias: 



e§-^Mt) = I{t) 



(27) 




FIG. 2: Time evolution of local electron density ni(t) in the 
central quantum dot in the system of three coupled quan- 
tum dots if energy level position in the middle quantum dot 
depends on time. Black line corresponds to the case when fre- 
quency of energy level ei (t) switching is lower than tunneling 
rates and grey line corresponds to the case when frequency is 
higher than tunneling rates. 



One can find for < t < to 



I{t) = eriijr 
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Mean tunneling current value can be found as: 
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If ri = l/2<o ^ 7res tunneling current mean value can 
be written as < I >— eVl since = I for such frequen- 
cies. This is the regime, when the device operates like 
a current standart: the current is directly proportional 
to the gate frequency. This regime has exponential accu- 
racy which is governed by the second and third terms in 
the square brackets in expression ([29| . Note, that even 
for very unsuitable case if |£2 — esl — 7 — ^ the pumping 
effect still remains, and the cuurrent is proportional to 
the frequency, though it's value is suppressed compared 
to the ideal relation < I >— eO.. 

For high frequencies of the gate voltage in the region 
7 3> O = I/2to 3> 7res tunneling current average value 
is almost independent on the frequency and equal to: 

< I >— e7res/4 (7res — ^f~)- With further frequency 
increase (f2 3> 7 3> 7res) mean current value decreases 
to < / 672^3/47. 

The non-stationary mean tunneling current value has 
non-monotonic dependence on the gate frequency with 
maximum at f2 ~ jres (FigE]). This effect can be used 
for frequency stabilization in nanoelectronics. 

The parameters of devices based on quantum dots 
depend on the size of quantum dots, tunneling trans- 
fer rates, energy levels positions and distances between 
them. Estimation of tunneling parameters for achievable 
setup gives us characteristic frequencies and currents for 
such devices: 
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FIG. 3: Frequency dependence of the mean tunneling current 
for the system of three coupled quantum dots (gate frequency 
Q, = l/2to). Different frequency scales are presented. 



T ~ 1 meV, 7 ~ 1 ^ 10 meV 

-fres = ^T^h - 0.1 ^ 1 meV ~ lO^o ^ 10" 1/sec 

where parameters T and 7 are determined by the 
widths and heights of the barriers. For such values of 
the tunneling rates we need to have quantum dots of 
tens of nanometers size, for which quantum size quanti- 
zation energies are not less than the tunneling width of 
levels. Such devices could operate at gigahertz and subgi- 
gahertz frequencies at nano and subnanoampere currents 
{InA ~ 6 • 10^ e/sec). 

The difference between tunneling rates in resonant and 
non resonant cases can be used also for creation of inverse 
occupation in quantum emitter based on the system of 
coupled quantum dots. An example of such device is 
shown in FigUl In the system of two coupled quantum 
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FIG. 4: Schematic diagram of two coupled quantum dots 
which operate as an emitter due to non stationary inverse 
occupation of levels. 

dots we switch the level e by applying gate voltage to the 
second quantum dot between two levels of the first quan- 
tum dot directly coupled to the lead. If the third non- 
resonant level exists in the same quantum dot between 
the lowest and the highest levels, then charge pumping 
from the initially filled state £3 to the highest level £1 
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creates inverse oecupation of the states £2 and Si. So 
some nanometer quantum generators can be fabricated 
on this principle. 

III. CONCLUSION 

We investigated electron pumping ability of a system of 
three tunnel coupled quantum dots attached to the leads. 
Periodical changing of energy level position in the middle 
quantum dot by gate voltage leads to nonzero tunneling 
current even if applied to the structure bias is equal to 
zero. 

Our calculations of the mean current are based on ac- 
curate analysis of relaxation processes in quantum dots in 
non-adiabatic regime. Exact equations allows to investi- 
gate various regimes of the device and estimate the mean 



current value and accuracy of it's operation as a current 
standard. For very small dots with pronounced size ef- 
fect a possibility of room temperature electron pumping 
is opened. We should like to stress that the ideas dis- 
cussed in the paper can be used for fabrication of a new 
type of electronic devices based on non-equilibrium non- 
stationary tunneling currents. As an example we pro- 
posed in the paper nanometer quantum emitter based on 
two coupled quantum dots. 
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